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5 We report the Fermi Large Area Telescope detection of 7-ray (>100 megaelec- 

e tronvolts) pulsations from pulsar J1823— 3021A in the globular cluster NGC 6624 

7 with high significance (~ 7 cr). Its 7-ray luminosity L 7 = (8.4 ± 1.6) x 10 34 ergs 

a per second, is the highest observed for any millisecond pulsar (MSP) to date, 

9 and it accounts for most of the cluster emission. The non-detection of the clus- 

10 ter in the off- pulse phase implies that its contains < 32 7-ray MSPs, not ~100 

11 as previously estimated. The 7-ray luminosity indicates that the unusually 

12 large rate of change of its period is caused by its intrinsic spin-down. This 

13 implies that J1823— 3021A has the largest magnetic field and is the youngest 

14 MSP ever detected, and that such anomalous objects might be forming at rates 

15 comparable to those of the more normal MSPs. 

1 6 Since its launch in 2008, the Large Area Telescope (LAT) on board the Fermi Gamma- 

17 ray Space Telescope (jT]) has detected whole populations of objects previously unseen in the 
is 7 -ray band. These include globular clusters (GCs), which are ancient spherical groups of ~ 
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10 5 stars held together by their mutual gravity. As a class, their 7 -ray spectra show evidence 
for an exponential cut-off at high energies ( 2 ][J|), a characteristic signature of magneto spheric 
pulsar emission. This is not surprising because radio surveys have shown that GCs contain 
large numbers of pulsars Q, neutron stars that emit radio and in some cases X-ray and 7 -ray 
pulsations. 

The first GC detected at 7 -ray energies was 47 Tucanae (5j), soon followed by Terzan 5 
and nine others (2j|3|). Even so, no individual pulsars in these clusters were firmly identified 
in 7 -rays Q. GCs are more distant than most 7 -ray pulsars observed in the Galactic disk (|<5j), 
thus most pulsars in them should be too faint to be detected individually. The Fermi LAT lacks 
the spatial resolution required to resolve the pulsars in GCs, which tend to congregate within 
the inner arcminute of the cluster. Hence, 7 -ray photons emitted by all pulsars in a given GC 
increase the photon background in the folded 7 -ray profiles of each individual pulsar in that 
cluster. 

One of the GCs detected at 7 -ray energies is NGC 6624 ([?), located at a distance d = 
8.4 ± 0.6 kpc from Earth (j9]) . With a radio flux density at 400 MHz of S400 = 16 mJy (10), 
J1823— 3021A is the brightest of the six pulsars known in the cluster. It has been regularly 
timed with the Jodrell Bank and Parkes radio telescopes since discovery, and with the Nan^ay 
radio telescope since the launch of the Fermi satellite. The resulting radio ephemeris (Table 
SI) describes the measured pulse times of arrival very well for the whole length of the Fermi 
mission, the root mean square of the timing residuals being 0 . 1 % of the pulsar rotational period. 

Thus we can confidently use it to assign a pulsar spin phase 0 to every 7 -ray (>0.1 GeV) 
photon arriving at the Fermi-FAT from the direction (within 0.8°) of the pulsar. We selected 
photons that occurred between 4 August 2008 and 4 October 2010 that pass the “Pass 6 diffuse” 
7 -ray selection cuts Q. The resulting pulsed 7 -ray signal (above 0.1 GeV, Fig. 1) is very robust, 


with an H-test value of 64 (11 ) corresponding to 6.8 a significance. The data are well modeled 
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by a power law with spectral index 1.4 ± 0.3 and an exponential cutoff at an energy of 1.3 ± 

0.6 GeV, typical of the values found for other 7-ray pulsars [see supporting online material 
(SOM)]. The two peaks are aligned, within uncertainties, with the two main radio components 
at spin phases (pi = 0.01 ± 0.01 and (p 2 = 0.64 ± 0.01 (Fig. 1). 

The pulsed flux above 0.1 GeV, averaged overtime, is F 7 = (1.1±0.1±0.2) x 10 -11 erg cm -2 s -1 , 
where the first errors are statistical and the second are systematic (SOM). The large distance 
of NGC 6624 implies that J1823— 3021 A is one of the most distant 7-ray pulsars detected 


([8]). This makes it the most luminous 7-ray MSP to date ( [72] ): Its total emitted power is 
L 1 = 47 rd 2 faF 7 = (8.4 ± 1.6 ± 1.5) x 10 34 (/h/0.9) ergs -1 . We obtained the statistical un- 
certainty by adding the uncertainties of d and F 7 in quadrature. The term fa is the power per 
unit surface across the whole sky divided by power per unit surface received at Earth’s location; 
detailed modeling of the 7 and radio light curves provides a best fit centered at 0.9, but with a 
possible range from 0.3 to 1.8 (SOM). 

The LAT image of the region around NGC 6624 during the on-pulse interval (0.60 < <p < 
0.67 and 0.90 < (p < 1.07) shows a bright and isolated 7-ray source that is consistent with the 
location of J1823— 3021A (Fig. [2]); in the off-pulse region (0.07 < <p < 0.60 and 0.67 < (p < 
0.90) no point sources in the energy band 0.1 - 100 GeV are detectable. Assuming a typical 
pulsar spectrum with a spectral index of 1.5 and a cut-off energy of 3 GeV, we derived, after 
scaling to the full pulse phase, a 95% confidence level upper limit on the point source energy 
flux of 5.5 x 10 -12 erg cm -2 s -1 . Thus, J1823— 3021A dominates the total 7-ray emission of 
the cluster. The combined emission of all other MSPs in the cluster plus any off-pulse emission 
from J1823— 3021 A is not detectable with present sensitivity. No other pulsars are detected in 
a pulsation search either. 

Under the assumption that the 7-ray emission originates from NGC 6624, ([?]) estimated the 
total number of MSPs to be TVmsp = 103 Assuming an average 7-ray luminosity for each 
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MSP (5 |2j), similar to the approximation made by Q, our off-pulse flux upper limit implies that 
/Vmsp < 32. This is consistent with the estimate TVmsp = 30 ± 15 derived from the correlation 
between 7 -ray luminosity and encounter rate (j2). Clearly, the MSP number estimate of (j+) is 
skewed by the presence of a single bright pulsar contributing disproportionately to its emission 


(13). The off-pulse emission limits can also be used to constrain alternative models for the 7 -ray 


emission from globular clusters, like those invoking inverse Compton (IC) radiation (14 15). 

The spin period of J1823— 3021A, 5.44 ms, is typical of MSPs. However, its rate of change 
-fobs = +3.38 x 10 - 18 s s 1 is one to two orders of magnitude larger than for other MSPs with 
the exception of J 1824— 245 2A, a pulsar in the GC M28 (16) that has a similarly large P D b s <Hzt- 
A possible explanation is that P 0 } )S is due mostly to the changing Doppler shift caused by the 
pulsar’s acceleration in the gravitational field of the cluster along the line of sight (a;): 

-fobs 


P I \ P c 


( 1 ) 


If the globular cluster has a reliable mass model, we could use it to estimate lower and upper 


limits for ai and estimate upper and lower limits for P (18). For NGC 6624 the collapsed nature 
of its core precludes the derivation of a reliable mass model. Furthermore, radio timing (Table 
SI) shows that J1823— 3021A is only 0"4 ± O'/I (a projected distance of 0.018 ± 0.004 pc) 
from the center of the cluster (79]), where the values of ai can be largest. For this reason, it 


has been suggested (10) that J1823— 3021A is a “normal” MSP (i.e., with small P ); its large 
P 0 bs being due to its acceleration in the cluster. This conclusion was apparently strengthened 


by the detection of a second derivative of the spin period P = —1.7 x 10 ^ s 1 (20). This 


could originate in a time variation of a/ resulting from interaction with a nearby object (21). If 
sustained it would reverse the sign of P 0 bs in ~ 6000 years; suggesting again that the large P 0 bs 
is not only due to dynamical effects, but is possibly a transient feature. 

However, the total observed 7 -ray emission L 7 must represent a fraction 7 < 1 of the 
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available rotational energy loss, E = r 2 PP/P 3 , where / is the pulsar’s moment of inertia. 
Although / depends on the unknown mass of the pulsar and the unknown equation of state for 
dense matter, the standard asumption / = 10 45 g cm 2 is a reasonable value for a 1.4-M 0 (mass 
of the Sun) neutron star. This implies P > 3.4 x 10~ 19 (/n/0.9)(//10 45 gcm 2 ) -1 ss -1 . Thus 
even an unrealistic y-ray efficiency rj = 1 would imply that P is already ~10% of P obs . If we 
assume instead P ~ P obs , thenP = 8.3 x 10 35 erg s _1 and rj = 0.1 x (/n/0.9)(//10 45 gcm 2 ) -1 . 
Comparison with the observed y-ray efficiencies of other MSPs ( 12 8) shows this to be a more 
reasonable range of values; r/ ~ 0.1 also represents the upper limit derived for the average effi- 
ciency of MSPs in 47 Tucanae (|5). Therefore, our y-ray detection of J1823— 3021 A indicates 
that it is unusually energetic and that most of P ob s is due to its intrinsic spin-down. The pulsar 
has other features that suggest it is indeed unusually energetic: Its alignment of radio and y-ray 
profiles has previously only been observed for the Crab pulsar (|22b and three fast, energetic 


MSPs: J 1939+2134 (the first MSP to be discovered), J1959+2048 (23) and J0034— 0534 (24). 
Like some of these energetic pulsars and PSR J1824— 2452A, J1823— 3021A emits giant radio 
pulses (25) and has a high 400 MHz radio luminosity of L 40 o — 1.1 Jy kpc 2 


), the third 

highest among known MSPs. However the correlation between E and radio luminosity is far 
from perfect given the uncertainties in the distance estimates, moment of inertia, beaming ef- 
fects and possibly intrinsic variations of the emission efficiencies. Finally, J1939+2134 also has 


a large P (26), which is thought to be caused by timing noise (TN), which scales roughly with 
P~ L1 P (27). In the case of J1823— 3021A, if P ~ P Q bs , then TN should be one order of mag- 
nitude larger than for J1939+2 134; instead its P is ~ 1.5 x 10 2 larger than that of J1939+2134. 
This is possible given the observed scatter around the TN scaling law. Thus TN might account 
for the P of J1823— 3021A, but this is far more likely if P ~ P obs . 

If P ~ P 0 bs, we can estimate the strength of its surface dipole magnetic field: P 0 = 3.2 x 


10 19 G\/PP(//10 45 gcm 2 )(P/10km) 3 ~ 4.3 x 10 9 G (28) [where Pis the neutron star (NS) 
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radius, generally assumed to be 10 km]. MSPs are thought to start as normal NSs with Bo ~ 
10 11-13 G which are then spun up by the accretion of matter and angular momentum from a 
companion star. This process is thought to decrease their magnetic field to Bo r\j 10 7 " 9 G; but 
the exact mechanism responsible for this is currently not well understood. Our value of B 0 
shows that for J1823— 3021A this decrease was not as pronounced as for other MSPs. 


As accretion spins up the NS, it eventually reaches an equilibrium spin period (29) given 


by: 


-Pin*. = 2.4ms 


6/ 7 / M ^ - 5/7 


(AG ) , 

Vio 9 g; \m ( 


f R \ 18/7 
\10 4 m / 


M 


- 3/7 


Mi 


(2) 


Edd 


where M is the NS mass, M is the accretion rate and M Edd is the maximum possible stable ac- 
cretion rate for a spherical configuration (known as the Eddington rate). Beyond this, the pres- 
sure of accretion-related radiation starts preventing further accretion. After accretion ceases, 
the newly formed radio MSP will have P init as its initial spin period. Assuming M = M Edd , 
M = 1.4M e andP = 10 km (as in our estimates of P 0 ), we obtain P init = 1.9ms(P o /10 9 G) 6//7 . 
For the value of B 0 calculated above, we get P init = 6.6 ms; that is, even if accretion had 
proceeded at the Eddington rate, the pulsar would not have been spun up to its present spin 


frequency. This is also the case for the other such “anomalous” MSP, J1824— 2452A (77); for 
all others we have P > P init . A possible explanation is that for these two objects M and / 
do not correspond to the assumptions above. If, for example, 77 = 0.15, M = 1.8 M e and 
I = 1.8 x 10 45 g cm 2 (30) we obtain B 0 = 3.6 x 10 9 G and P init = 4.7 ms. A second possibil- 
ity, suggested by eq. [2} is super-Eddington accretion (more precisely, M > 1.6 iif E dd); this can 
happen for non-spherical mass accretion. A third possibility is that the value of B 0 was smaller 
during accretion (resulting in a smaller Pinit), and that B 0 has increased since then. This has 


been observed for some normal pulsars (31 ); however there is no evidence of such behavior for 
any other MSPs. 
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In any case, the conclusion that P ~ P Q b s implies a characteristic age r c = P/ (2 P) = 25 
million years. This is likely an over-estimate of the true age of the pulsar, particularly given 
that P init is likely to be similar to P. Thus J1823— 3021 A is likely to be the youngest MSP ever 
detected; only J 1824— 245 2A might have a comparable age. Because of their large P s both 
objects will be observable as MSPs for a time that is ~ 10 2 shorter than the ~ 100 “normal” 
radio-bright MSPs known in GCs. Statistically, this suggests that, at least in GCs, anomalous 
high 73-field MSPs like J1823— 3021A and J1824— 2452A are forming at rates comparable to 
those of the more “normal”, radio-bright MSPs. 
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Pulse Phase 


Figure 1: Phase-aligned radio and 7 -ray profiles for J1823— 3021A. (Bottom) Nangay 1.4 GHz 
radio profile. (Top and middle) 7 -ray profiles obtained with the Fermi-LAT in different energy 
bands. The dark histogram is for events with E > 3.0GeV. The 7 -ray background for the 
0.1 GeV light curve was estimated from a surrounding ring, and it is indicated by the dashed 
horizontal line in the top panel. The highlighted area there shows the on-pulse region selection. 
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Figure 2: Fermi LAT 7-ray count map above 100 MeV for J1823— 3021A during the on-pulse 
(Left) and off-pulse (Right) regions, as defined in Fig. 1. The 6° by 6° region is centered on 
the pulsar position (cross). The map was adaptively smoothed by imposing minimum signal- 
to-noises ratios of 13 and 16 for the on- and off-pulse regions, respectively. The large circle 
indicates the tidal radius of NGC 6624. The small circle shows the 99% confidence region for 
the location of the 7-ray source. 
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Supporting Online Material 

Observations and data analysis 

Radio Timing analysis 


With the express purpose of supporting the Fermi mission (33), J1 823— 3021 A is observed 


approximately 3 times per month with the 76-m Lovell telescope (34). using a 64 MHz band 
centered at 1404 MHz connected to an analog filterbank. Since mid-2009 observations have also 
been performed using a digital filterbank backend with 1024 x 0.5 MHz channels of which ap- 
proximately 250 MHz is used. Highly precise timing measurements are also conducted with the 


Nan<jay radio telescope (35). These have included regular observations of J1823— 3021A since 
mid-2006. Approximately every two months, the pulsar is observed for 1 hour at 1.4 GHz. 
A 128 MHz bandwidth is coherently dedispersed using powerful GPUs (Graphics Processing 
Units). A total of 104 pulse times of arrival (TOAs) were obtained from the two telescopes 
between mid-2006 and mid-2010. The Tempo2 timing package (36) was used to build the 
timing solution, which includes the pulsar rotation frequency and its derivatives, the dispersion 
measure, and the pulsar position. The post-fit residuals are characterized by a weighted rms of 
7.3 /is. The resulting parameters are summarized in Table SI. No trends are noticeable in the 
post- fit residuals. 


Fermi LAT data analysis 

We have observed J1823— 3021 A with the Large Area Telescope aboard Fermi from 2008 Au- 
gust 4, when the satellite began scanning-mode operations, to the end of the validity range of the 
pulsar radio ephemeris (2010 October 14). The data analysis presented in this paper has been 
performed using the LAT Science Tools package 09-21-00 and the P6 V3 Diffuse instrument 
response functions (IRFs). Events tagged “Pass 6 diffuse” having the highest probability of 
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being 7 -ray photons (37) and coming from zenith angles < 100° (to reject atmospheric 7 -rays 
from the Earth’s limb) were used. Additionally, a rotational phase was assigned to each selected 


LAT event using the radio ephemeris as an input to the Fermi plugin (55) distributed with the 
Tempo 2 pulsar timing software. 

Using the pyLikelihood likelihood fitting tool with the NewMinuit optimizer, we per- 
formed a binned spectral analysis to determine the energy flux and the spectral shape of the 
source. Events in the range 0. 1-100 GeV were extracted from a 20° x 20° square region of in- 
terest (ROI) centered on the pulsar position. To reduce the effect of the Earth’s atmospheric 
emission, the time intervals when the Earth was appreciably within the field of view (specifi- 
cally, when the center of the field of view was more than 52° from the zenith) were excluded 
from this analysis. The Galactic diffuse emission was modeled using the gllJem_v02 map cube, 
while the extragalactic emission and residual instrument backgrounds were modeled jointly by 
the isotropic component isotropic _iem_v02. These two models are available from the Fermi 
Science Support CenteiQ In addition, all the sources found in an internal catalog based on 18 
months of data (similar to (| 59[)) above the background with significances > 5a and within 20 ° 
from the pulsar were included in the model. Sources were modeled with a power law spectrum, 
except for pulsars for which a power law with an exponential cut-off was used (4C |. Sources 
more than 5° from the pulsar were assigned fixed spectra taken from the source catalog. Spectral 
parameters for sources within 5° of the pulsar were left free for the analysis. 

The pulsar location at the core of NGC 6624 is located just outside the 99% statistical 


error contour of the 7 -ray source 1FGL J1823.4— 3009 (39), based on an analysis of the first 


1 1 months of the LAT survey data. Therefore, (41 ) did not establish an association between 


1FGL J1823.4— 3009 and the globular cluster. (42), using a larger dataset, showed that the 
7 -ray position lies within the (20.55') tidal radius of the cluster. To check the association, 


1 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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we first reevaluated the position of the 7-ray source 1FGL J1823.4— 3009, using the on-pulse 
(0.60 < <p < 0.67 and 0.90 < <p < 0.07) segment of the pulsar rotational phase (to improve 
the signal-to-noise ratio), to be (a 2 ooo = 275.87°, 5 2 ooo = —30.29°) with a 99% confidence 
error radius of 0.09°. This places the pulsar radio position just inside the 99% error contour 
of the 7-ray source. We also investigated the off-pulse window and detected a ~ 4 a point 
source which coincides with the position of the radio source NVSS J 182324— 30031 1 (a 20 oo — 
275.8515 °, 5 2 ooo — —30.053278 °). The signal, only observed during the first year of the Fermi 
mission, is located 0.1 1 ° away from the position of 1FGL J 1823.4—3009 and 0.31 ° away from 
NGC 6624. It is likely that 1FGLJ 1823.4-3009, located between NVSS J182324-30031 1 and 
NGC 6624, includes contributions from both of these sources. Using the second year of data, 
when the nearby faint source is off, we localize the 7-ray source corresponding to the pulsar to 
(^2000 = 275.93 °, 5 2 ooo = —30.34°) with a 68% error counter radius of 0.07 °. This position is 
consistent with the radio pulsar position. We then fitted the spectrum of the pulsar at the radio 
pulsar position using a power law with an exponential cut-off. Figure S 1 shows both the fit 
between 0. 1 and 30 GeV (solid lines) and the spectral points derived from likelihood fits to each 
individual energy band in which it was assumed the pulsar had a power-law spectrum. 


Light Curve Modeling 


There have been two major contenders for modeling the high-energy (HE) radiation (roughly 
100 MeV to 10 GeV range) from pulsars, those which assume that the observed 7-rays are emit- 


ted near the stellar surface (43 44) above the magnetic polar cap and those which assume the 
7-rays originate primarily in the outer magnetosphere near the light cylinder (the distance from 
the rotation axis Q at which the co-rotation equals the speed of light). Both classes of models 
assume that the HE 7-rays are curvature radiation from highly-relativistic electrons/positrons in 
the radiation-reaction regime. The two most common outer-magnetospheric emission models 
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are the outer gap (OG; (45, 46] p 7])) and the slot gap (48) models. For our purposes we took the 


two-pole caustic (TPC; (49)) model to be a geometric realization of the slot gap. TPC and OG 
models both assume that the emitting electrons are accelerated up to high altitudes in narrow 
gaps along the last-open field lines. The OG model only allows acceleration above the null- 
charge surface (NCS; where Q ■ B = 0) whereas in the TPC model electrons are accelerated 
from the stellar surface. The HE pulse profiles in these outer-magnetospheric models are the re- 
sult of the accumulation of photons in narrow phase bands due to a combination of three effects: 
aberration (change of photon direction due to the high corotation velocity), time-of-flight de- 
lays (photons produced at higher altitudes will reach an observer earlier than those coming from 
lower altitudes), as well as the magnetic field line curvature (photons are assumed to be created 
tangential to the local magnetic field line in the corotating frame, and their direction is therefore 


very sensitive to the magnetic field geometry). This is referred to as caustic emission (50). 

TPC and OG models are generally used in conjunction with a low-altitude radio cone beam 


geometry (e.g., (51 52)). Due to the difference in altitude of the radio and y-ray emission, there 
will be a phase lag between the radio and y-ray profiles. Polar-cap (e.g., (44)) y-ray emission 
models do predict much smaller phase lags but, due to the large open field line region of MSPs, 
they cannot produce the narrow peaks observed in the y-ray light curve of PSR J1823— 3021A. 
The phase-alignment of PSR J1823— 3021A’s radio and y-ray light curves argues for overlap- 
ping y-ray and radio emission regions. To reproduce the phase-aligned light curves we used 
altitude-limited versions of the TPC and OG models (alTPC and alOG, respectively) which 


were first introduced to model the light curves of the MSP PSR J0034— 0534 (53j). These are 
very similar to the standard TPC and OG models, except that the minimum and maximum radii 
of the radio emission region as well as the maximum radius of the y-ray emission region are free 
parameters (the minimum y-ray emission radius being set by the standard models). Therefore, 
both radio and y-ray photons originate in a TPC or OG-like structure, with a significant amount 
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of overlap between the two emitting regions leading to phase-aligned profiles. This implies 
that the radio emission is also caustic in nature, supported by polarimetric observations which 


find 0% linear polarization for PSR J1823— 3021A (54). Conversely, these models provide a 
framework to constrain the respective radio and 7 -ray emission geometries when comparing the 
model light curves to the data. 

We have simulated 7 -ray and radio light curves using alTPC and alOG models with a spin 
period P = 1.5 ms, steps of 1° in magnetic inclination angle (a) and viewing angle (0, 0.05 
in accelerating emission layer width (w, normalized to the opening angle of the polar cap), 
and 0.10 (in units of f?LC = cP/(2tt)) for the emission altitudes. The spin period used in 
the simulation is less than that of PSR J1823— 3021A (5.44 ms) but this quantity enters the 
simulation mainly through the size of the polar cap. Using models with a shorter period will, at 
most, overestimate any predicted off -pulse region. We have developed a Markov chain Monte 
Carlo (MCMC) maximum likelihood technique to jointly fit the 7 -ray and radio profiles and 


pick the best-fit paramters (55). We fit the > 500 MeV 7 -ray light curve in 60 bins and rebinned 
the radio profile to 60 bins. For the 7 -ray models the minimum emission altitudes (R^J were 
specified as described previously while the maximum emission altitudes were allowed to be free 
under the constraint R)) l0:J . > 0.7 Rlc- The radio emission altitudes are unconstrained save that 


Rjj )r(r > R^ in . The best-fit parameters for both models are given in Table S2; the likelihood 
does not prefer one model over another. The alTPC model has best-fit gap widths of 0.0. This 
is unphysical and should be taken to indicate that the true gap widths are between 0.0 and 


0.05. Following (52) we can estimate the beaming factor (fa) for both models using Eq. 4 
of (57) , see Table S2. Presently we are unable to provide reliable uncertainty estimates for our 
model predictions and, thus, can not propagate any uncertainty on fa into the uncertainty on 
L 7 . However, while PSR J1823— 3021A stands out in some respects, the shape of the observed 


HE light curve is very typical of known 7 -ray MSPs. The best-fit geometries of (55) to these 
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448 MSPs yield values of from approximately 0.3 to 1.8, with mean of 0.81 and rms of 0.36. We 

449 therefore expect the /q value for J1823— 3021A to be similar to the geometries which, in these 

450 models, produce “typical” 7-ray light curves. The high and low tails of this distribution suggest 

451 that the 7-ray efficiency could reasonably be anywhere from 3 to 20% but neither extreme 

452 affects the conclusion of the main text that most of the observed P is intrinsic to the pulsar. 
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Timing parameters 


Right Ascension, a (J2000) 18 h 23 m 40:48466(4) 

Declination, <5 (J2000) -30° 21' 39"988(4) 

Solar System Ephemeris DE 405 

Reference time scale TDB 

Reference time (MJD) 54939 

Spin Frequency, v (Hz) 183.823389814514(7) 

First derivative of u, v (10“ 15 Hz s -1 ) .. —114.1351(4) 

Second derivative of v, i) (10~ 25 Hz s -2 ) 5.8(1) 

Dispersion Measure, DM (cm~ 3 pc) 86.864(9) 

Validity Range (MJD) 53773.35-55483.67 

RMS Timing Residuals (/is) 7.3 


Table SI. Timing parameters for J1823— 3021A. The center of the globular cluster is located at 

a = 18 h 23 m 40!51 ± (k008, 5 = -30°21 / 39'/7 ± O'.'!. 



Fig. SI. y-ray spectral energy distribution of J1823— 3021A obtained with the Fermi Farge 
Area Telescope. The solid black line shows the maximum likelihood fit to a power law with 
exponential cut-off. The dashed lines are ±1 a uncertainties on the fit parameters. Plotted points 
are from likelihood fits to individual energy bands with > 3<r detection above background for 
two degrees of freedom, otherwise a 95% confidence level upper limit arrow is shown. The 
errors are statistical only. 
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Fig. S2. Observed and best-fit 7 -ray (top) and radio (bottom) light curves for 
PSR J1823— 3021A using the alTPC (pink) and alOG (green) models described in the text. 
The dashed, horizontal lines in both panels correspond to the estimated background levels. The 
7 -ray background was estimated using an annular ring centered on the radio position with inner 
and outer radii of 1° and 2°, respectively. The radio background was estimated by fitting the 
region between 0.1 and 0.6 in phase to a constant value. The parameters of the best-fit light 
curves are given in Table S2. 
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